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Abstract 



We consider the interface problem between two incompressible and inviscid fluids with constant densities in the presence of surface 
tension. Following the geometric approach of [14, 15] we show that solutions to this problem converge to solutions of the free-boundary 
Euler equations in vacuum as one of the densities goes to zero. 

1 Introduction 



r Z~ 1 1.1 Description of the problem and main results 
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We consider the interface problem between two incompressible and inviscid fluids that occupy domains Qf and fl^ in R" (n > 2) at 
time t. We assume fijj" is compact and M™ = Q,f U fl^ U St where St := dflf. We let v±, p± and p± > denote respectively the 
velocity, the pressure and the constant density of the fluid occupying the region ttf. We assume the presence of surface tension on 
the interface which is argued on physical basis to be proportional to the mean curvature k + of the hypersurface St- The equations of 
motion are given by 1 

' p(v t + v ■ Vv) = -Vp S t 



V ■ v = x g R™ \ St 

v(0,x) = v°{x) x G M™ \ S 

with corresponding boundary conditions for the interface evolution and pressure's jump given by 

d t + v± ■ V is tangent to \J t S t C R n+1 
p+(t, x) - p-(t, x) = e 2 K + (t, x) , x e S t ■ 



(E) 



(BC) 



We are interested in analyzing the asymptotic behavior of solutions of the above equations when p_ — > 0. Our result is convergence 
to the solution (v+, S^°) of the free-boundary problem for Euler equations 



p+(d t v + + v + ■ Vu 4 
V • v + = 

t>+(0, x) = v+(x) 



x G f2 t °° 



x G 



(Eo) 



with corresponding boundary conditions 



where k°° denotes the mean curvature of S? 



d t + v+ ■ V is tangent to (J t S%° C M" +1 

P + (i,a;) = 6 2 K°°0,a;) , a; G 5 t °° 

:= dVlf. More precisely we will show the following 



(BC ) 



'Here we are introducing the notation / = /+Xq+ + f—Xn— ^ or an y /± defined on f2 t . 
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Theorem 1.1. Let an initial hypersurface 1 So G H l+l and an initial velocity field vo G H l (R n \So) be given for some I > § + 2. 
Consider any sequence of local in time solutions of (E)-(BC) 

S?£C([0,T],H l+1 ) , v m €C([0,T},H l (Q™)) 

corresponding to values of the density p m = p+Xn+ + Z 9 - Xn~ w ^ tn P™ ~ > ^ as m ~ > 00 • Then (v™, S™) converge^ on a small time 
interval to the solution 

St = dSlf e C(H l ' +1 ) , v+ G C(H l '(Q%°)) far any I' < I 
of(Eo)-(BCa). Convergence is in the space (S t ,v) G L°°(H l -i) x L°°(H l ')for any I' < I. 

Free-boundary problems for Euler equations have been extensively studied in recent years following the breakthrough of Wu in [ 1 8, 1 9] 
where local well-posedness in Sobolev spaces is proved in 2 and 3 dimensions for the irrotational gravity water wave problem. Many 
works have dealt with the water wave problem also in the general non-zero curl case, see for instance [13, 8, 14, 9]. 
For the irrotational vortex sheet problem with surface tension Ambrose [1] and more recently Ambrose and Masmoudi [2] proved 
well-posedness respectively in 2 and 3 dimensions. Cheng, Coutand and Shkoller [6] proved well-posedness in 3-d for the full 
problem with rotation and well-posedness is also obtained (in any dimension) by Shatah and Zeng [16] for (E)-(BC) and other realted 
fluid surface problems [16, sec. 6]. 

In absence of surface tension the vortex sheet problem for the free-boundary motion of two fluids is ill-posed due to the Kelvin- 
Helmotz instability as shown in [11]. Beale, Hou and Lowengrub [4] showed how the surface tension regularizes the linearized 
problem. In the next section we will show how the Kelvin-Helmotz instability is very apparent from the infinite-dimensional geometric 
arguments presented by Shatah and Zeng in [15]. 

We recall that also the free-boundary problem for Euler equations in vacuum (Eo)-(BCo) with e = is known to be ill-posed due to 
Rayleigh-Taylor instability, see [10], which occurs if one does not assume the sign condition 

-V N+P+ (x,t) >a>0 . (RT) 

In [14] it is shown how also the Rayleigh-Taylor instability is a natural consequence of a geometric calculation and is related to the 
sign of an operator appearing in the linearization of the Euler flow. Motivated by this we are going to show 

Proposition 1.2. Let T be the space of all admissible Lagrangian maps for the interface problem (E)-(BC) defined in (1.7) and let 

r* := {$ : Qq — > W 1 volume-preserving homeomorphisms} 

be the corresponding space for the water wave problem (Eo)-(BCo). Consider a point u G T and tangent vectors Vi G T U T for 
i = 1 . . . 4, where T U T is endowed with the L 2 (p m dx) metric. If we denote' 1 ' TZ and TZ* the curvature tensors ofT and V* respectively, 
then 

(K(u)(v u v 2 )v 3 ,v 4 ) L 2 {pmdy) "=5° (TZ*(u + ){v 1+ ,v 2+ )v 3+ ,v 4+ ) L2{p+dy) (1.1) 

In view of the geometric frame work described below and the linearized equation (1.17), proposition 1.2 can be considered as a first 
step in showing that solutions of (E)-(BC) converge to solutions of (Eo)-(BCo) with e = when e,p- —> at the same time 5 . 

Our paper is organized as follows. The geometry of T is presented in section 1.2 and an explanation of the geometric intuition behind 
the Kelvin-Helmotz and Raileigh-Taylor instabilities is given in 1.2.3. Of course we refer to [14, 15] for full details about this general 
geometric approach. In section 2 we state theorems on energy estimates which are independent of p_ . Proofs are performed in section 
3 and appendix B. Section 4 is devoted to showing strong convergence of solutions as stated in theorem 1.1. The proof of proposition 
1.2 is then performed in section 5. 

During the writing of this manuscript it was brought to the attention of the author that Cheng, Coutand and Shkoller [7] had proved an 
analogous result to the one stated in theorem 1.1. 

1.2 The geometric approach to Euler equations 

It is well-known that the interface problem between two fluids has a variational formulation on a subspace of the space of volume- 
preserving homeomorphisms. For the water wave problem this was observed for the first time by Arnold in his seminal paper [3], 
where he pointed out that Euler equations for the motion of an inviscid and incompressible fluid can be viewed as the geodesic flow 
on the infinite-dimensional manifold of volume-preserving diffeomorphisms. This point of view has been adopted by several authors 
in works such as [17, 5, 12] and more recently by Shatah and Zeng in [14, 15, 16]. 

2 The regularity of hypersurfaces in R™ is intended in the sense of local coordinates: an hypersurface is H s for s > 2 if it can be locally represented as the graph 
of H s -functions. 

Convergence is achieved by reducing the problem to the fixed initial domain Qq using Lagrangian coordinate maps. See section 4 for details. 
4 Covariant differentiation on T U T (and on T U T*) is defined in section 1.2.2. 

5 We believe that some condition of the form p_ = 0(e a ) for some a > should be needed in this case. 
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1.2.1 Lagrangian formulation 

The surface tension parameter e will be henceforth set to be one. Multiplying (E) by v, integrating over K"\5t, using the boundary 
condition (BC) and the variation of surface area formula, we obtain the conserved energy 6 

E = Eo(S t ,v)= [ ^-dx+ f dS=: [ ^-dx + S(S t ). (1.2) 

For y G Qq we define u±(t, y) to be the Lagrangian coordinate map associated to the velocity field v±, i.e the solution of the ODE 

^=v±(t,x) , x(0,y)=y VyellJ; (1.3) 

for any vector field w on R n \St we define its material derivative by 

D t w := Wt+v ■ Viu = (w o u) t o . 
In [15, sec. 2] the authors derive from (E)-(BC) the following equation for the physical pressure: 

-Ap = ptr(Dv 2 ) 

P±\ St = M- 1 {-^M TKT -2V v T_ v TV^-Tl + (vl,vi)-Tl-(vI,vI) (1.4) 

-Vw.A^trpt) 2 ) - V N _AZ 1 tr(Dv 2 )} 

where n ± denotes the second fundamental form of the hypersurface S t (with respect to the outward unit normal vector N± relative to 
the domain Clf) and Af is given by 

p+ p- 

with N± denoting the Dirichlet-to-Neumann operator on the domain Q t . 

From (1.3) we see that in Lagrangian coordinates Euler equations assume the form 

pu*t = -Vpou w(0)=idn o (1.6) 

with p determined by (1.4). 

Since v is divergence free, u± are volume-preserving maps on M™\5 , . Moreover u + (t, So) = U-(t, So) even if the restriction to So 
of u + and h_ do not coincide in general. This leads to the definition of the space T of admissible Lagrangian maps for the interface 
problem: 

T = |<E> = $+Xn+ + ®-Xn~ '■ *± : ~ * *±(^o ) is volume-preserving homeo. and d$±(fij) = $±(90q ) j . 
Denoting S(&) = J$(s ) ^ we can rewr i te me ener gy (L2) in Lagrangian coordinates as 

E {u,u t )= f ^^dy + S(u) 
where p = pou. The conservation of the above energy suggests that (E)-(BC) has a Lagrangian action 



// 

J Jr 



P\ut\ 



R"\S ^ 



I (u)= / / ! - L ^ j -dydt - / S(u)dt. (1.7) 



6 Notice that the conserved energy does not control the L 2 norm of V- in the asymptotic regime p_ — » 0. 
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1.2.2 The geometry of T 

To derive the Euler-Lagrange equations associated to the action I one has to look at the geometry of V considered as a submanifold of 
L 2 (pdy) and identify its tangent and normal spaces. It is easy to see that the tangent space of T at the point $ is given by divergence- 
free vector fields with matching normal components in Eulerian coordinates 7 

T^T = jw) : R"\5 -> K" : V • w = and w+ + w_| <&(So) =0, where w = w o . 
while the normal space is 

(r$r) ± = {-(v^)oc f > : kcso) = p-V>-|* (So) = : V> s } • (i-^) 

A critical path •) of J satisfies 

©t«t + 5'(u) = (1.9) 

where 5'(u) denotes the tangential gradient of S(u) and D t is the covariant derivative on T (along u(t)). In order to verify that the 
Lagrangian map associated to a solution of (E)-(BC) is indeed a critical path of (1.7) we need to compute S' and D t . 
Computing <D t an d H u {t)'- Given a path u(t, •) G T denote v — u t and St — u(t, So). For any w(t, •) G T U T we must have 

w t = V t w + II u{t) (w,v) (1.10) 

where II u (t) {w, v) G (T u ( t )r) denotes the second fundamental form on T u ( t )I\ From (1.8) there exists a unique scalar function p ViW 
defined on M. n \S t such that 

II u{t) (w,v) = -Vp v , w oue {T U ^T) ± (1.11) 

In [15] it is shown thatp ViW is given by 8 

—Ap v ^ w = tr (DvDw) 

Pt, w \ St = ^, tu --^- 1 {v„T^T«;|+V u) T_ tu Ti;|-n + («J,^) (1.13) 

-n_(«I,ioI) - Vjv+A" 1 b(DvDw) - X7 N _ A: 1 tr (DvDw)} . 
Then in Eulerian coordinates we can write 

D t w := (i> t w) o u^ 1 = D t w + Vp v ^ w . (1.14) 
Computing S'(u): For any w G T U T the formula for the variation of surface area gives 

{S'(u),w) L2(Rn ^ SQpdy) = K+W^dS 

J St 

and it is not hard to verify that the unique representation in Eulerian coordinates of S'(u) as a functional acting on T U T is 

S\u) = Vp fe with p± = -L-UtM^M^ . (1.15) 

P-P+ 

From (1.4), (1.13) and (1.15) we obtain the identity p = p(pk + p v ,v)- Therefore, taking w — u t , we see from (1.14) and (1.15) that a 
solution of (1.9) equivalently satisfies 

T> t v + Wpv, v + Vp K = (1.16) 

which is exactly (1.6) in Eulerian coordinates. 



7 We follow the convention used in [15] where the Lagrangian description of any vector field X : &(Qq) —* is denoted by X = X o $. 
8 Let us point out that in the water wave problem with just one fluid in vacuum we have II*^(w, v) = — VpJ w ou€ (T u ^r*) ± with 

-&P*v,w = tr (DvDw) 

(1.12) 

Pv,w\gn t — 0- 
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1.2.3 Linearized equation and instability 

The Lagrangian formulation discussed above provides a convenient setting to study the linearization of the problem. Considering 
variations around the solution u t of (1.9) and taking a co variant derivative with respect to the variation parameter, we obtain the 
following linearization for w{t, ■) 6 T U ^T: 

i 2 w + %.(u)(u t ,w)u t + i 2 S{u)w = (1.17) 

where denotes the curvature tensor of the manifold T and D 2 S(u) is the projection on T U ^T of the second variation of the surface 
area. Both of these linear operators acting on T U T play a central role in the understanding of the problem and in the definition of 
high-order energies based on their leading order terms. In [14] an explicit but rather complicated formula is given for D 2 S(u); in 
[14, 15] its leading order term j? is singled out and turns out to be given 9 in Eulerian coordinates by 

*(u)M=V/+xn++V/_xn- with f± = -^—^Af^M^-AsJwi ; (1.18) 

P+P- 

it is easy to see that j? is a third-order 10 self-adjoint and positive semi-definite operator with m(u)(w, w) = | Vw^ \\^(s t )- Further 
computations performed in [15, pp 859 - 860], show that the leading-order term Ha{u){v) of the unbounded sectional curvature 
operator %. (u) (v,-)v is given in Eulerian coordinates by 

*,(«)(«;) = V/+xn+ + V/_xn- with /± = -^—^N^N^^ ^N^V ■ (wi{vl - vl)) . 

P+P- 

Noticing that % (u) is a second-order negative semidefinite differential operator we immediately see that the linearized Euler equations 
would be ill-posed if there had been no surface tension generating the operator a . This is the so-called Kelvin-Helmotz instability for 
the two fluids interface problem. 

We conclude this section recalling that the same geometric setting described above applies to the problem of Euler equations in vacuum. 
The same Lagrangian approach is of course available and the linearized equation is still given by (1.17). Computations performed in 
[14, sec 2.2] show how the leading order term of the differential operators involved in the linearization are given by ^* (u) and a * (u) 
satisfying 

ni(v,w) = H(q {u) + bounded operators 

D 2 S(u) — * (u) + second-order differential operators 

and 

^(u) = [ -V NP * vv \Vw^\ 2 p+dS , a*(u) = [ IVw^p+dS. (1.19) 
Js t Js t 

Since also in this case A*{u) is generated by the presence of surface-tension, we see that (1.17) is ill -posed in absence of surface 
tension if the Raileigh-Taylor sign condition (RT) is not assumed. 

2 Theorems on Energy Estimates 

Definition 2.1. Let Ao = Ao(So, I — \, S, L) for some I > ^ + 1, L > and < S <C 1 be the collection of all hypersurfaces S such 
that a diffeomorphism F : So — > S CK" exists with 

\ F - id ^„'-i ( s 0) <5 
and satisfying a uniform bound on the mean curvature \k\ ,_b - < L. 

In [15] the geometric considerations exposed in section 1.2 led the authors to define the following energy for (E)-(BC) 

9 Both in the one fluid case and the interface problem the leading order term of i> 2 S(u) has the same form but its Hilbert space representation does not coincide due 
to the different orthogonal splitting of L 2 in T$r and (T$T) X in the two settings. We refer to [15, pp. 857-858] for the details of the derivation of m. 
'"Assuming St is smooth enough. 
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Definition 2.2. Consider domains Clf with Clf compact and interface St G H l+l . Letv(t,-) G H l (K™\ St) be any divergence-free 
vector field with v + + v_ = 0, define the energy 

If L 2 If L _i 2 2 

E(S t ,v(t,-)) = - \*3v\pdx+- |*s *Vp K \ pdx+\u\ H ,- 1(R n^ St ) (2A ^ 

w/zere is z7ze cmW o/u, f/iaf is or- = 9^ — c^V. 

Proposition 2.3. Lef / > | + 1, f/zen for S t G A vvz'f/z S t G i? (+1 we have 

\k\ 2 h1 - 1{s) <C (1 + E) , \v\ 2 Hl(R ^ s) <Co(l + E + E ) 2 

where C is some positive constant depending only on A and the initial data ( in particular it is independent of p_). 

The above proposition is the equivalent of [15, proposition 4.1]. The proof of bounds which are independent of the density p- just 
requires some small modification of the argument given in [15]. See section 3. 

Theorem 2.4 (Energy Estimates for (E) and (BC), [15]). Let I > \ + 1 and a solution to (E)-(BC) be given by 

S t G H l+1 and v G C° (H l (R n \ S t )) , 

then there exists L > and a positive time t* independent of p- and depending only on \v(0, OlijirR^xSt)' an d L, such that 
St G Ao and |«|^fi-i(5 t ) < L for allO <t < t*. Moreover the following energy estimate holds for < t < t*: 

E(S t ,v(t,-)) < 3£7(S ,«(0,-)) + Ci+ [ P(E ,E(S t ,,v(t',-)))dt' (2.2) 

Jo 

where P is a polynomial with positive coefficients determined only by A and the constant C\ depends only on \vq\ ; _ a and 
A . 

The proof of theorem 2.4 is essentially the same as in [15] and is postponed to the appendix. 
Corollary 2.5. Consider a sequence of solutions 

S™ G C° (H l+1 ) , v m G C° (iJ i (M"\S* t m )) 
solving locally in time the Euler system (E)-(BC) for values of the density p™ — * 0. If we denote 

E m (t) :=E(S?,v m (t,-)) 

with E given by (2.1), then there exists a positive time and a constant C depending only on the set A , | w oIh ! (K"\S ) an d 
\ v o\ H i-% {Rn ^ So) such that 

sup E m (t) < 2E(0) + 2Ci , V m G N. (2.3) 

te[o,tj] 

The above corollary gives as a consequence weak convergence of solutions of the vortex sheet problem to solutions of the one fluid 
problem in vacuum. Standard compactness arguments are going to give the strong convergence stated in theorem 1.1. See section 4 
for details. 

For completeness we state here a theorem, proved in [16], based on the above energy estimates and concerning existence of solutions: 

Theorem 2.6 (Well-posedness for (E)-(BC), [16]). Given an initial surface So G H t+1 and initial velocity v Q G H l (W l \ So) with 
I > § + 1, the free interface problem (E)-(BC) has a solution in the space 

S t G C° (H l+1 ) , veC° (H l (R n \ St)) 

for t in some small interval [0, T] independent of the density p_. Moreover, if I > 3 the problem is locally well-posed, i.e. the solution 
is unique and depends continuously on the initial data. 
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3 Proof of Proposition 2.3 

Using the definition of A in (1.18) we can explicitly write the terms appearing in the energy (2.1) as in (B.2), (B.3) and (B.4) with k/2 
replaced by 1/3. From the properties of A/± and A/" -1 in lemma A. 2 it follows that there exists a constant C independent of p_ such 
that 

\K\ 2 Hl - liSt) <C(l + E) , \v±f Hl - i{St) <C(l + E). 
To estimate v we proceed in three simple steps: 

1 ) Estimates on the Lagrangian coordinate map: Consider the Lagrangian map U- associated to V-. From lemma A.l we get 

id| H . (n -) ^ Cl J Q l v -(*.-)lH.(n t -)l t *-(*.-)lff.(n 5 -) dfl 

for any < s < I where C\ > only depends on n and I. Now, let ^bea sufficiently large number to be specified later depending 
only on the _ff'-norm of the initial velocity, define 

t :=sup{i : \v(t',-)\ Hl{Rn ^ Sti) <» Vi'G[0,t]}. (3.1) 

Since v is assumed to be continuous in time with values in H, to > 0. The previous inequality and an easy bootstrap argument (or 
Gronwall's inequality) show that there exists a positive time t\ and a constant C*2 depending only on Z, n,fi and A such that 

|«_(*, ■) - id n - < C 2 t < ^ , Vte[0,t*] (3.2) 

for t* :— min{t , 1/(2C2)}. This shows that w_ is an iJ z -diffeomorphism so that •) is a well-defined volume preserving 

map for x G fi^~ and for the same range of times we have 

K^-rV(n-) <2 , v o<*</-i. 

2) Decomposition of vector fields and control of the lower norm: As it is well-known (and explained in detail in [14, Appendix B]) 
any divergence-free vector field v : R n \S t — > R" obeying the condition w^: + = can be decomposed into two divergence-free 
components, the rotational part v r responsible for the interior motion and an irrotational or gradient component Vi r = \7g responsible 
for the motion of the boundary S t - More precisely g is the solution of the elliptic Neumann problem 

| Aj = i6l\S, 
1 VAr ± .g± =v^ x £ S t 

and v r := v — v ir . It is observed in [16] that the invariance of Euler equations under the action of the group of volume preserving 
diffeomorphisms leads, via Noether's theorem, to a family of conserved quantities which determine the rotational part of the velocity 

v r (t, •) = P r (S u (Du-'fviO, u~\t, •))) (3.3) 

where P r (S tl X) denotes the projection of X : R n \S t — > R™ onto its rotational (gradient-free) part. Therefore we can estimate 

MV)l L 2(n-) < M L 2 (0t -) + K(Mli,2(n-) < \(D U - 1 )*v(o, u - 1 (t,-)\ L2{n - ) + \v-(t,-)\ H i {St) 

< KDurX^Jv^^y^+CEi <C(1 + E$) 
with C depending only on the initial data and Ao. 

3) Control of \v\ Hi : To conclude we use the fact 11 that any divergence-free vector field can be controlled by its curl and normal 
component: 

\v\% ( nt) < ^a + l^+l^-f^^O^^-Mo^ + Kl^-^so + l^^)) <C(1 + E + E f 
where the constant C depends only on the initial data and the set A □ 

1 'A more general statement is 

\ v ±\Hi(n±) < + l K +l H i-| (s) ) (\ divv ±\H'-Hn ± ) + \ cmlv ±\H'-Hn ± ) + \ v ± - VA_1 divu±| ffI _ i ^ + \v±\ L 2 (n± ^j 
where the constant C depends only on An. An essential proof of this can be found in [14, proposition 4.3]. 
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4 Proof Theorem 1.1 

In this section we are going to use the uniform bounds provided by corollary 2.5 combined with the non-linear Eulerian frame work 
introduced in [ 14] to obtain the strong convergence of solutions stated in theorem 1.1. 

4.1 Convergence of Lagrangian maps and velocity fields 

As a first step we need to estimate the physical pressure. 

Lemma 4.1. Let v G H 1 and S t = dfl t G H l+1 with I > ^ + 2 be a given solution of (E)-(BC). Then the pressure p, determined by 
(1.4), satisfies 

^H'-hnf) ^ C (\ V +\H'-Hnt) + \ K +\w- H s t )+P-\^-i(^sJ N ^-Hs t) ) (4-1) 

and for p_«l 

\P-\ H '-h ( n-) ^ °P- {\ V -\ 2 H'-Hn-) + l«-lff.-i( St) + V-xsJ^-MS*)) (4-2) 
for some constant C depending only on the set of hypersurfaces Ao 

Proof Write p± — A^ 1 Ap± + H±p± and use lemma A.2 to get 

< Cp±\v±\ 2 H ,- 1(nt) + C^\K±\ H ,- 1(St) + Cp- (\N±\ Hl _ 2(St) \v\ 2 Hl _ i{Qt) 

Proposition 4.2. There exists a sequence {nik}, a time t** depending only on the initial data and an H l -diffeomorphism u + G 
C? ([0,t**];H l (n+)) withd t u+ G C° ([0,t**];i? z '(f7+)) suchthat 

limu™*=u+ in C t °(([0,t**];iJ'(fi+))) (4.3) 
Urn <9 t u™ fc = d t u+ in C° ([0, t**]; ^'(0+)) (4.4) 

for any I' < I. Moreover if we define 

:=u+(i,fi ) (4.5) 
?nen f/zere exists u+ G i°° (H l {n™)) n L°° smc/z f/zaf 

lim u™* o = w+ o w+ in C t ° f [0, ***]; H 1 ' (Q^)) (4.6) 

for any I' < I andp + G L°° (H l ~i (ftf)) sucn mat 

lim p r l h o u" lfc =p + ou + weak-star in L°° ( [0, i**]; #'"5 (4.7) 
We will still denote these subsequences by the index m. 

Proof Let us denote X(H S ) = A([0, t**], H S (Q^)) for X = L°° or C t ° and C any positive constant depending only the initial data 
and the set A . Combining proposition 2.3 and corollary 2.5 we see that 

\ vm \L°°(H'(ci+- m )) ^ ^0(1 + Em) < C 

for any t < Therefore, arguing as in the proof of proposition 2.3, we can find a positive time t** < t$ depending only on A and 
the initial data, such that for any < t < t** 

- id 0+ | . < Ct** < i. (4.8) 
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This show that each map w™ is an TJ'-diffeomorphism onto its image and is uniformly bounded in L°°{H l ) by a constant depending 
only on the initial data and the set A . Then, up to extraction of a subsequence, there exist u + G L°°(H l ) such that u™ — > u + 
weak-star in L°°(H l ). Lemma A.l and (4.8) imply 

lft<l a . ( nf)<Wl a .(# n l<t( W f,<o. 

Again by standard compactness we have, up to extraction, <9 t w™ = v™ o u™ — > <9 t u + =: 5+ weak-star in L°°(H l ). Since u + G 
W 1 ' 00 ^'), we get u+ £ C t °(iJ') and 12 u™ -> u+ in C?(JT'). 

Passing to the limit in (4.8) we see that u + is also an ff'-diffeomorphism. Thus we can define v + by v + o u + =: v + = d t u + . From 
Euler equations we have d t (u+ o u™ ) = — Vp™ o u™ so that lemma 4.1, lemma A. 1 and corollary 2.3 together with (A. 1 0) imply 

In particular this gives continuity of v+ o u™ = dtu+ with values in i/' _1 (f2[j~). It also implies the existence of a subsequence (still 
denoted by the index m) such that d t (v™ o u™) — > weak-star in L°°(H l ~i). Since u" 1 oii™^ij + in the sense of distributions, 
V+ = d t v+. Therefore 13 v+ = v+ o u+ e C^H 1 ' 1 ) and 

v™ou™^v + ou + in C?(([0,t**];^' _1 M)) • 

As i>+ o u" 1 is uniformly bounded in L°°(H l ), by interpolating the Sobolev norms we can improve the above convergence obtaining 

(4.6) and the equivalent (4.4). 

Finally, since p™ o u™ is uniformly bounded in L°°(H l ~ 5 ), up to extraction, we have p™ o u™ —> p+ weak-star in L°°(H l ~ 2 ) and 

(4.7) follows just by defining p + =: p + o u^ 1 □ 

4.2 Verification of (BC ) 

Using convergence of the Lagrangian maps u™ associated to u" 1 established in (4.3), we defined in (4.5) the "limit domain" where 
the evolution of the limit solution is going to take place. From (4.3) and trace estimates we obtain w+| s — > M +ls ^ C®{H l ~^ (So)) 
so that 

u+(t,Sb) = 5u+(t,0+) =: £ C?(ff'-') for * G [0, t**] . 
Proposition 4.3. 77;e moving boundary condition in (BCo) /ioWi /or f/ie set of hypersurfaces with v + defined by (4.6). 

Proof From the definition of Lagrangian maps, (4.4) and (4.6) we have 

d t u + (t,y)=v + (t,u + (t,y)) V (t,y) G [0,t**] x 0+ . 

As u + (t, So) = S^ 00 for any i G [0, t**], we have that (t, u + (t, •)) is a curve on the space-time boundary U t S^°; therefore 

d t + d t u+ ■ V = d t + v+ o u+ ■ V is tangent to (J C M™ +1 . 

t 

The fact that u + is a diffeomorphism from So to for any t G [0, i**] gives the claim □ 

Lemma 4.4. Lef N™(t, •) fee f/ze outward unit normal and k™ (t, •) the mean curvature of S™. Denote by N°°(t, x) and n°°(t, x) 
respectively the unit normal and the mean curvature o/5 t °° at the point x. Then for any V < I 

N" 1 o u ™ — > N°° o u+ in C t °(i? r {S )) and K ™o«^ K M o w+ in C?(H'' _1 (S )) • (4-9) 

/n particular -1(500) is uniformly bounded which implies 14 S^° G if +1 as stated in theorem 1.1. 

12 The standard argument is the following. Consider an arbitrary subsequence of {«+}; the boundedness of {9t«™} implies through the Ascoli-Arzela theorem the 
existence of a sub-subsequence converging in C^(H l ) to a limit which must be it+ (the weak * limit of the original sequence {u™}). Therefore is the uniform 
limit of {u™}. 

13 We use the fact that / g L 2 {H^ ) and ft € L 2 (H S ^) imply / € C(H^+ S ^/ 2 ). 

14 This can be proved using local coordinates and estimates for quasi-linear elliptic equations. Another proof can be found in [14, proposition A.2]. 
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Proof Since k°° (t, x) (X, Y) = tr (Y ■ V X N°° (t, x) ) for any X, Y G T X S^°, it is enough to prove the first statement in (4.9). 
We use similar arguments to those in the proof of proposition 4.2. By lemma A.l, (A. 10) and (2.3) we obtain uniform bounds on 
Af™ou™ in L°°(H l ); therefore there exists N + G such that, up to extraction of a subsequence, N™ o u™ — > N + =: A°°ou + 

weak-star in L°°(H l ). Identity (A. 6) and estimate (A. 10) combined with the uniform energy bounds on k™ show that 



^ ) <CK\^ ) \N-\^ ) <C (4.10) 



with some C uniform in Ao and m. This in particular implies that N™ o u" 1 belongs to C(H (So)) and that, up to further extraction, 

d t {N™ o u") -► dt(A°° o u+ ) weak-star in L^H^iSo). As a consequence, A°° o u+ G C^tf'" 1 ^)) and N™ o u ™ -^A oc ou+ 
in C t °iJ r (So) for any I' < I. 

To show that A°°(t, •) is the outward unit normal N?°(t, •) to the hypersurface let r m £ T X S™ be an arbitrary tangent vector. 
Since u™ is a diffeomorphism from So to S^ 1 , there exists a unique tangent vector t € T'yS'o such that r m = rfu™(t, j/)ro, where 
du™(t, y) denotes the differential of as a map from So to S™ acting on T y So for y = (w™)~ (t, a;). Then for any t G [0, t**] 

(N? (t,u™(t,y)) , du™(t,y)r ) =0 

Letting m go to infinity using (4.3) we obtain 

(A°° {t,u+{t,y)),du + {t,y)T ) = 

Since r m , and consequently to, was arbitrarily chosen this implies that A°°(t, x)A.T x S£° for x = u+(t, y)\ by the strong convergence 
established above ^4°° is unitary and therefore coincides with N°° (t, x) □ 

Proposition 4.5. The boundary condition (BCo) for the pressure is satisfied by the limit solution. 

Proof For the sequence of solutions (u m ,S™) condition (BC) holds for every m G N. As (E) is also satisfied for every m, the 
boundary condition for the physical pressure is the one given in (1.4) (where of course every quantity has to be indexed by m). 
Therefore (p™ — k™) o u™ = p™ ° w+ on So and we can use lemma A.l, (4.2) and trace-estimates to obtain 



IW-0 <l ff «- 1(So )<^ ro | H --i . . • 



Since the expression in parentheses above is uniformly bounded by the energies, letting m — > oo and using (4.9) we get 



>k°°o U+ in r; //'' '(S,,) (4.11) 



P+ ° u+ 

for any /' < I. Using (4.7) we conclude that p + (t, x) = K°°(t, x) for any t G [0, t**] and x G S t °° □ 

4.3 Verification of (E ) 

We first need the following estimate: 
Lemma 4.6. Letp + be given by (1.4) then 



^\ L oo (H ,- a(n p m)) <C. (4.12) 



for some C uniform in m. 



Proof In what follows we suppress the use of the index m and let a < b denote a < Cb for some constant C independent of /?_ . 
Writing p + = H+p+ + A -1 tr (Dv + ) 2 we have 

T> tP+ = U t H +P+ + D t A 1 tr (Dv+ f = H+B tP+ + A ^ tr (Dv+ f +R:= (I) + (II) + R (4.13) 

where the remainder is given by the sum of the two commutators 

R = R 1 +R 2 := [T> t ,H+] p+ + [D t) A" 1 ] tr {Dv+f . (4.14) 
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We show that every term is bounded in H l ~ 2 or better by the quantities \v\ Hl , \p\ H i- 1 > Njji-i an d \N\ H i which are already known 

to be bounded uniformly in time by the energies independently of p_ . 

Estimate of (I): This is the highest order term in (4.13). Denoting P := Afp+\ s we have 

(I) = H+DtM^P = n+Af^BtP + n+R^P with R 3 := [A/" _1 ,D t ] . 
Observe that R 3 = Af~ x [Af, D t ] AT' 1 so that (A.3), (A.5) and (A. 13) give 

l W + i?3p l ff '-i ( o+) £ \W^^] p \w- H s t )^P-\ v \w(n t) \P\H^(s t) 
< \v\ H ' { n t ) (l«lff«-i ( s t ) + Mff'-i^J^-MSt)) • 

Using again (A.3) and (A.5) we obtain 

\H + Af-i-D t P\ Hl _ 2{nt) < Cp.\D t P\ Hl _ i{St) . (4.15) 
Now D t P contains four different terms to be estimated. The term involving the mean curvature is estimated by (A. 13) and (A. 7): 

(\v\ Hl{nt) \K+\ Hl _s {St) + \ V+ \ Hl{nt) \N + \ Hl _ i{St) + l«+l ff .-i (St) KI H «-i (n +)) • 



< ^ 

~ P- 



Notice that the presence of p_ in the denominator in this last estimate is compensated by the factor in (4.15) so that the bounds 
remain uniform. For the terms involving tr(Dw) 2 we use (A.6), (A.12) and the identities D t V/ = VD t / — (Dv)*Vf and 15 
D t tr (Dv) 2 = -2 tr [(Dv) 3 - 2p+D 2 p ■ Dv] to estimate 



iDtV^A^tr^i) 2 ^,.^ < \n t N ± \ Hl . HSt) \vf Hl _s {nt) + iJVil^^jlDtVA^trpVi) 2 ^, 



'(fit) 

) 



< |^V±| H! -2 (St) b|^_i (at) + |JV±|ffi- a(lSt) (l«±l^ ( ^)l w ±l^-Mo t ) + l«±l!r'-i(n t ) + IP±l ff ^ ( o t ) |w lff , - 1 (o t ) 

Analogously, using (A. 8) the terms T) t Il±(v±, v±) and D t (v±Vv±) can be bounded uniformly in H l ~^(S t ) and H l ~ 3 (S t ) respec- 
tively. 

Estimate of (II): By the same formula used above to express D t tr (Dv+f we get 

lA-^tr^) 2 ^.^ < IWVi^+A-H^ 

£ (KI^- 1 (o+) + b+l^-l ( a+)l t; +l^-Mn+)) • 
Estimate of R: Commutators R 1 and i? 2 are estimated directly by (A. 11) and (A.12): 

|[D t ,W+]p+| Jf! _ 2(a + ) < b+l^(o+)l^+l^-i(o+) ' [D t ,A- 1 ]tr(^+) 2 + £ l«+lHi (n +) 
where as usual the constant C is independent of p_ □ 



15 This identity follows from H t Dv = DD t v — (Dv) 2 together with Euler equations pD t u = — Vp. 
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Proposition 4.7. Let v + and u + be given as in proposition 4.2 then 



l^ou™)— ►|(t, + o« + ) in c?(H''-t(n+)). 

and v + satisfies Euler equations (E ). 

Proof (4.6) and the uniform bounds on p 7 ? establish the above convergence weak-star in L°°{H l ~i {£1$)). Since d?(v™ o u™) = 
D t Vp!p o u™ = VD ( p™ o u™ - (D^J^Vpf o the bound given in (4.12) implies <9 t 2 (v™ o u™) e L 00 ^" 3 ^)) and the 
desired strong convergence follows through the usual arguments. 

From (4.7) and (4.1 1) we know that p" 1 o u" 1 — > p + o u + strongly in C^H 1 ~^ and therefore Vp™ o u+ = V(p+ o u™)(V)t™) 1 — ► 
Vp + o m + in C i °i? i_2 (fi^). Since Euler equations in Lagrangian coordinates are <9 t (w™ o = — Vp™ o u™ we can take the limit 
in L 00 (i? i_2 (ilQ )) obtaining that v+ satisfies Euler equations in Lagrangian coordinates too, i.e. 

j t v + (t,u + {t,y)) = -Vp + {t,u + (t,y)) V (t, y) e [0, t**] x fi+ . 

Finally from (4.3) and (4.6) we have V{v™ o u™) -> V(w+ o u+) in C* t °(ff'- 2 (f7+)) so that 

= V • o u ™ = tr(Vw™ o u ™) = tr (v« o «;)(V<) _1 ) ? ™ tr (v(v+ o u+ )(Vu+) _1 ) = V ■ v+ o u+ 

which implies V • v + = pointwise in for any t G [0, t**] □ 
The proof of theorem 1 . 1 is completed ■ 

5 Proof of Proposition 1.2 

Let r be the infinite-dimensional manifold (1.7) and TZ its curvature tensor induced by the co variant differentiation defined in section 
1.2. Consider a map u(t) : Qq — > O t m T. Let 7?™ denote the sectional curvature of T at the point u as an operator acting on T U T 
endowed with the L 2 (p m dy) metric and depending on some v e T U T (and of course on u). We assume v and the hypersurfaces St to 
be sufficiently smooth and single out the leading order term of TZ m analyzing its behavior as m goes to infinity (or equivalently as the 
density p_ vanishes). 

In view of the geometrical frame work discussed in section 1.2, and in particular in 1.2.3, TZ m can be considered as a measurement of 
the instability occurring in the linearized Euler equations in case surface tension were not present. 

Let w be any vector in £ T U T. We assume that w is uniformly bounded in H l (M. n \St) for some large enough I and compute the 
sectional curvature in the direction of v and w. Using a well-known formula from Riemannian geometry together with (1.11) we have 

TZ m = CR(u)(v,w)v,w) L2(pmdx) = {II u (v,v), II u {w,w)) L 2 (pmdx) - \\IIu(v,w)\\ 2 L2(pmdx) 

= / yp v ,vVpw,w P m dx \Vp v ,w\ 2 P m dx . 

Again we suppress the use of the index m. Using the divergence theorem the first integral can be written as 

/ Vp v . v Vp w . w pdx = Pvv^n+Pww + ^N-Pww) p v . v Ap w . w pdx 

Jl»sS ( J S t jR"\S t 

= f P S vv {-2V w t twI +Tl + (wl,wl)+Tl_(wI,wI)\ dS+ f p v , v tr(Dwfpdx 

J S t L ' JR"\S t 

having used Vn + Pw. w + ^N-Pw.w = ^Pv.v + ^ 7 w+A^ 1 Ap+ UJ + Vjv_ Al 1 Ap~ w and (1.13) with v = w. Since tr(Dw) 2 = 
diW k dkW l — di(w k dkW l ) we can use twice again the divergence theorem obtaining 

/ ^Pv.v^Pw.wpdx = / D 2 p ViV (w,w) pdx + p% v |-2V ffi T_ ffi Tw{ + II + (tt;+, «;+) +U-(wI,wI) 

jR"\S t jR"\S t JSt L 

+ X7 w+W+ -N + + X7 w _w_-N^} dS- [ p+w^\7 w+P + v +p_w^\7 w _p- v dS. (5.1) 

Js t 
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To estimate the terms containing p^ w , which is the inverse image through J\f of a mean zero function on St, we use lemma A. 2. For 
any / G L 1 (S' t ), (A.5) yields 



f plJdS <C\^ V J \f\ Li{St) <Cp-W^J H .-, {st) \f\L 1 

J St 



whenevers > 2—^, with C uniform in S t <E A . Since \NPv, w \ Hs -i( S -j is uniformly bounded for smooth enough and bounded v and 
w, we can easily estimate several terms in (5.1): 



< Cp-\w\ 3 , 



J St J St J St 

where C is some uniform constant depending on v and the mean curvature of S t . These bounds imply 

lim / Vp v . v Vp w . w pdx = / D 2 p v+iV+ (w + ,w + ) p + dx - lim / p + w+V w+ p+ v + p-W±V w _p~ v dS . 
p-^° jR"^St Jn+ P-^°JSt 

Next we look at the contribution of \\II u (v, w)\\ 2 , use the decomposition f± = H±,f + A^ 1 A/ applied to p VyW and observe that 
VH±±VA± 1 A to obtain 

/ \Vp v , w \ 2 pdx = / p^ w Afp^ w + / |VA; 1 tr(Z?wL>w)| 2 / 9 + dx+ / IVAI 1 tr {DvDw)f p_ . 

In [15] it is shown how the leading order term of the sectional curvature comes from the contribution of the surface integral in the 
above expression and is a second order negative semi-definite operator. But since Afpy w is independent of p_, by the same argument 
performed above the boundary integral vanishes as /?_ — > 0. Therefore 

lim K m = / D 2 p vv (w,w)p+ - \V AT 1 tr ( DvDw)f p+dx 
- lim / p + w+V w+ p++p-W±V w _p-dS. 

P-^ a JSt 

By splitting w into normal and tangential components on the boundary the surface integral in (5.2) is 

/ P+wi V w+ p+„ + p-w±V w _p~ dS = / p + \w^f\7 N+ p+ v + p-\w^f\7 N _p- V 
J s t J s t 

+ / P+w^Vl p+ + p-w^:VZ_p~ dS . 
J s t 

Writing p±pf v = H± (Py v ) — p±A^} tr {Dv) 2 , by the usual estimate foxp^ v the right-hand side of (5.3) gives the contribution 

/ K f^+P'L + \ w - f-^-PL + P+ \< | 2 V W+ A; 1 tr (Dv) 2 + p_ | 2 V N _ AI 1 tr (Dv) 2 dS 
J s t 

Pz ^° I p+V N+ A^tv(Dv) 2 \wi 
J s t 



(5.2) 

(5.3) 
(5.4) 



dS. 



(5.5) 



Since p+p^ v = P-P v v = Py v on S t and we are considering only tangential derivatives, the contribution of the term in (5.4) is 

/ W^\7 w T +w T P l v dS < C\w\ 2 Hl pn^ St) \p$ iV \ H , l{St) < Cp-\ W \ 2 Hl{Rn ^ St) "^"0. 



Gathering (5.2), (5.3), (5.5) and (5.6) we get 



lim 

m— >oo 



(K m (u)(v,w)v,w) LHpmdx) + [ P+ V N+ A-^(Dv + ) 2 \wi\ 2 dS 

J St 



< C\w\ 



(5.6) 



(5.7) 
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This shows that the leading order term of the sectional curvature of T in the limit p_ — > is given by the self-adjoint operator lZo(v) 
acting on T u+ T represented in Lagrangian coordinates by 

n (v+)= (-P+VH+ (v N+ A-hr(Dv + ) 2 (- \ du+{no) f)) ou + 

and satisfying 

(TZo(v + )w + , w+) L2[p+dy) = - [ Vjv + A- 1 tr(^+)Vi| 2 p+^- (5.8) 

From (1.12) we see that A -1 tr (Dv + ) 2 is exactly p* + v+ for the water wave problem so that (5.8) is equivalent to the first integral 
in (1.19). Therefore we have shown that as p_ — > the Kelvin-Helmotz instability for the vortex-sheet problem becomes the 
Raileigh-Taylor instability, i.e. the leading order term of the sectional curvature of V is not definite in general and has a positive sign 
only provided that the normal gradient of the physical pressure (in absence of surface tension) is negative. We conclude with two 
observations: 

(i) If we do not restrict our attention exclusively to the highest order term of the sectional curvature operator, the above calculations 
show that 

Jfa (Tt m (u)(v,w)v, w) L2{pmdy) = - f p+VN+A-hriDv+flw^dS 

- / D 2 A- 1 \.t(Dv + ) 2 {w + ,w + ) P+ -\V A- 1 tr (Dv + Dw + )\ 2 p + dx. 

Since the second fundamental form of T* in the water wave problem is given by Vp* = VA _1 tr (DvDw) the above limit is 
exactly 

/ VPv + ,v + Vpl +iW+ P+dx- / \Vp* v+ , v+ \ 2 P+dx= (H i '(u + ){v + , w+)v+, w+) L 2, dy) . 

(ii) From a standard argument we conclude that the full curvature tensor of T converges to the curvature tensor of T* in the sense 
stated in (1.1) and this completes the proof of proposition 1.2 ■ 

A Supporting material for proofs 

In this appendix we gather some technical results needed in the proofs presented and in the proof of theorem 2.4 in appendix B. 

Lemma A.l. Let D t (resp. Si) be domains (resp. hypersurfaces) in W 1 for i = 0, 1. Let r\ : D — > -Di (resp. r\ : So — > Si) be an 
H l -diffeomorphism for I > ^ + 1 (resp. I > I1 ^-)with | (det D-q)^ 1 | L oo( Dl ) < a (resp. | (det D , q)~ 1 \ Laa ^ Si ^ < a). Then the operator 
Tr) ■ / — > / or] is a bounded operator from H s (Di) to H S (D ) (resp. from H S (S\) to H s (So)) for any s E [0,1] and satisfies 

\f°ri\ H °(D ) < Co (uoMffi^) ( A -!) 
for some constant Co depending on a, s, I and the domains Di (resp. the hypersurfaces Si). 

Proof The case s = follows immediately from the hypotheses. Assume by induction that (A.l) holds for any integer s such that 
0<s<fc-l<Z-l. We prove the statement for s = k. Write D k (f orj) = D k ~ 1 (Df o r\ Drj) = J^jZo D j (Df o 7 1 )D k -3i 1 . Let 
r > 2 be the integer such that ^ — 1 < Z — r < observe that D tr q e L°° for i < r — 1 while it is not uniformly bounded in general 
for i > r since H l ~ l does not embed in L°°. According to this we split 

fc— 1 k—r k — 1 

^D j (Df or])D k - j r] = Y, Di ( D f or])D k - j Ti+ D j (D.f o V )D k - : >f] =: Sj + S 2 . 

j—0 j—0 j—k — r+1 

In S 2 a H derivatives on rj can be taken in L°° and estimated through Sobolev's embedding: 

fc-i 

|S 2 | L2(Do) < I^W ^)lL 2 (Do)l^"^lLoc( Do )<^/^l^- 1(Do )W ff «Po) 
j— k— r+1 

< C \ D f\ H ^-HD 1 )\ r l\m(D )\v\m(D ) = C \f\ H * (Di) I 7 ? \H l (D ) • 
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The contribution of Si is estimated using Holder's inequality and Sobolev's embeddings. Since I > ^ + 1 and k — j > r > 1, we can 
choose 2 < p, q < oo such that 

1 1 _ 1 1 1 l-k + j 1 1 k-l-j 
p q 2 ' q 2 n p 2 n 

Using Holder's inequality and the embeddings #'- fe +J c ff^-f )« c £,9, ftk-x-i c #(3-f)™ c £ P we get 

fe— r 

I S iIl 2 (Do) < ^^(Wor?)!^^)!^"^!^^) ^ C f |£^'(D/o^)| Hib _ 1 _ ; , (Do) |^| Hl(Do) <q/| jrfc(Di) |^|* rI(Do) 

with C depending on a, k and the domains Do, Di. Therefore we proved 

\DHf°v)\ LHDo) <C\f\ Hk{Di) \v\ k Hl{Do) 

From the inductive hypothesis the same inequality holds for lower order derivatives terms D l ( f o rj), for < i < k — 1, replacing fc 
with i. Up to further increasing the value of C depending on a and the constants in Sobolev's embeddings, we can sum this inequalities 
to obtain (A.l) for s = k. The case for non-integer s follows by interpolation □ 

Lemma A.l (About differential operators on A [14]). Let A -1 and H denote respectively the inverse Laplacian with Dirichlet 
boundary condition and the Harmonic extension 16 operators on a domain 0. There exists a uniform constant C > such that for 
every domain ft with <9£1 = S e A (see definition 2.1) the following is true: 

\fU HH on) < C \f\ H * + i m - Vs>0 < A - 2 > 

|A-V-(n),^ +1( n)) + \ n \ HH ^ m) .H^ m ^ C > V fl €(0,f-1]. (A3) 

As a consequence the Dirichlet-to-Neumann operator satisfies 17 

^LiH'+hen^-honv + ^ V.-i<«,),iz- + i<«,)) ^ C ■ V « € [0, / - 1] , (A.4) 

where H s denotes zero-mean H s -functions. In particular if Af is the operator defined by (1.5) then \-^\ L ^ H s+^^ dn ^ H s ~i(dn)) — 

C{P- +P+)/(p-P+) and 

W-X^-iw^hw) * 2 °P- > V a € [0,/ - 1] W p_ < ^ . (A.5) 

Proof The proof of (A. 2), (A. 3), (A.4) and more detailed analysis of operators acting on <9£1 (and in particular of the Dirichlet-to- 
Neumann operator) can be found in [14, A.2]. To obtain (A.5) write j\f as 

\P+ ) P- P- 

Estimate (A.4) implies that for p_ < p+/(2C 2 ) the linear operator B maps H s (dQ) to itself with norm less or equal than \. Thus 
/ + B is invertible and A/"" 1 = p-Ml 1 J2T=o B j s ° that 

OO 

3=0 

Lemma A.3 (Geometric Formulae [14]). Let S be an hypersurface in R™ moved by the normal component of a vector field v. Let 
N, k and II denote respectively its unit normal, mean curvature and second fundamental form. Then the following identities hold true: 

B t N = -[(Dv)*-N] T (A.6) 
B tK = -A St v- L -v- L |n| 2 + V„TK (A.7) 
BJU(t) = -V T ([{DvfN) T ) -n((V r «) T ) (A.8) 

-A s n = -2? 2 K + (|n| 2 /-/«n)n (A.9) 



l6 F = A^ 1 / satisfies AF = / in Q and F = on dCl. G = Tig satisfies AG = in tt and G = g on dfl. 

17 In view of (A.3) jV can be defined for any / g H s (dCl), s > \ in the weak form {<p, jV (/)) = f n VHipVHf. 
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where V denotes the covariant derivative on S and As := trV 2 . Furthermore there exists a uniform constant C such that for any 
SeA a 

|n| H . (s) + |JV| ff . +1(s) <c(i + |K| ff . (s) ) vz-|<*<i-i. (A.io) 

Lemma A.4 (Commutator Estimates [14]). There exists a uniform constant C such that for any dfl = S G Ao the following 
estimates hold: 

l[ D '' W V-W(S)) ^ vl< S <l (All) 

l^'^L^-w^)) < CMhw V2-;< s <Z (A.12) 

\\Pt,Mo]\ L(H . (S ),H-ns)) ^ C \ v \m { n) Vl<,s</-i (A.13) 

|[D t ,A s ]| L(ff3(s) ^_ 2(s)) < C\v\ Hl(n) -l <s <l-\. (A.14) 



B Proof of Theorem 2.4 



This section is devoted to the proof of Theorem 2.4 and consists essentially of material contained in [15, sec. 4.3, 4.4]. The only 
difference is that we claim and show independence of the energy estimates on the densities of the two fluids. Therefore, even though 
the proof is extremely similar to the one performed in [15], we present it here for the reader's convenience. 

B.l Estimates on the Lagrangian coordinate map 

We use the same notation in the original proof of theorem 2.4 letting I := |fc. Working on the compact domain Clf and arguing as in 
the proof of proposition 2.3 (see section 3) we obtain the existence of a positive time ti and a constant C\, only depending on k, n and 
(i as in (3.1) such that 

!«+(*.-)- id n + \ H iH (nt) <Cit V te [0,min{i ,*i}]. 
This implies the estimate on the mean curvature 18 

\K+(t,-)\ Hik _5 {St) <C 2 t+\K + (0,-)\ H 3 k _5 {So) V te [0,min{t 0) ti}]. (B.l) 

where the constant C2 is only determined by fj, and the set A . We conclude that there exists a time t 2 again determined only by \i and 
the set A such that 

S t GA , V te [0,mm{t ,t 2 }}. 

B.l Evolution of the Energy 

The energy defined in (2.1) is made of three terms. The first two involve the operator a defined in (1.18) and are used to control the 
irrotational part of the velocity and the mean curvature (hence the regularity of the evolving domain St); the third part involves the 
vorticity ui and is used to control the rotational part of v. More explicitly E = E\ + E 2 + |w|^3 fc _ 1 where, using (1.18), the first two 
terms are given by 

Ei := \ ! \^v\pdx = \( vl{-A St M) k -\-A St )vidS (B.2) 
E 2 := If \^--^ VK \ p dx = \\ K + Af(—As t J\f) k ~ 1 k + dS (B.3) 

1 JR"xS t 1 J S t 

where 

N= (^-Af^Af- 1 (J-JS) ■ (B.4) 
It is clear from lemma A. 2 that TV" is a first-order self-adjoint operator whose norm and inverse's norm do not depend on p_. 



18 This can be checked using the local coordinates constructed in [14, appendix A]. 
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Proposition B.l. There exists a polynomial Q(t) — Q ( \v(t, -)| 3 fc Mt, -)\ 3 k _ 1 ) with positive coefficients depending 

on the set Ag and independent of the density p- such that 



clt 



(E - E a ) 



< Q 



(B.5) 



where the extra energy term E ex , due to the Kelvin-Helmotz instability, is given by 

E ex = T( P l s I V v T K+ .tf(-A St tf) k -\ v TK + dS- P ~ 

2{P++P-)Js t + + 2 {P+ + P-)Js t 



\7 vI K + -M(~A St M) k 2 V vI K+dS. (B.6) 



Proof Throughout the proof we denote by Q any generic polynomial satisfying the properties in the statement. 

Evolution of E\. This is the hardest term to deal with and is the one where the extra energy term E ex appears. We are going to show 



d_ 
~dt 



dS 



<Q. 



(E 1 -E a )+ f «i(-A 5t A0V 

From definition (B.4) and (A. 13) we have 

[T> t ,Af] \ L(H s (SthH s-i (St)) < C f |«| ff4 » (fJt) ■ 2 _ - _ 2 - 2 
Therefore, using also (A. 14), we can commute D t+ with the operators appearing in (B.2) to get 



(B.7) 



n , V I < s < Ik - - 



<Q. 



~ j s vi{-A St M) k -\-A St )vidS- j s vi(-A St Af) k -\-A St )T> t+ vidS 

Using (A.6) to express D t+ N + , (1.16), (1.13) and (1.15) together with (B.4) we have 

T>t+v+ = D t+ v+-N + + v+T> t+ N + = -V N+ p+ v -V N+ pt -V v tv+-N+ 

= -—N+p s v v + V N+A+ 1 K{Dvf - Nk+ - V vI vl + U{vl,vl) 
p+ ' + 

From lemma A. 2 and trace-estimates \Vn + A^ 1 tr (Dv) 2 \ H 3 k _ i ^ < Q so that this term is lower order and 

jEi- j s vl(-A St Sr) k -\-A St )(-±Af + pl v -tfK + -V v T + vl + n(vl,vl^ dS 
Equation (1.13) forpf gives 

-^f+p S v , w = --^/v+AT 1 [2V v T_ vl vi -U + (vl,vl) -U_(vl,vl) -Vat+A; 1 ti(Dvf -V^AI 1 tr^) 2 } . 



< Q. 



(B.8) 



P± 



Since A/+7V 1 is an operator of order zero the terms Vtv± A^. 1 tr (Dv) 2 can be treated as before. From Lemma 4.6 in [15] (to which 
we refer for the proof) 



(-A St )^-A4 
this and the definition (1.5) of TV yield 

Together with (A. 10) this gives 



-±AT +P ? V — ^- 
P+ ■ P+ + P- 



< C (l + \n(t, -)| 3 k -i ( st)) V - --k < s <-k-- 



N+N- 1 



P+P- 



P+ + P- 



<Q. 



L(ffi fc -i(5 t )^3 fc -7(5 t )) 

(2V„t _ vI v$- n+ (vj ,vl)-n-(vl,vlj) 



H? k -i(s t ) 



< Q 
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so that (B.8) becomes 



d F 
Jt El 



vl(-^sM~\-^s t ) 



p- 



l( P-~ Pa 



P+ + P- 
2P- 



-U_(vl,vl)+ P+ - '« T 



P++ P 



-n+K,<) 



\p+ + p- p+ + p- 



dS 



< Q. 



We now claim that the last two terms in the above integral are lower order. To see this, consider flows < &±(t, •) on D,f generated by 
H+v± and apply (A. 13) and (A. 14) to D T to move outside the tangential derivatives V^t : 



/ «i(-A St A7) fe - 1 (-A 5t )V^ 

J St 



's 

p- - p+ 



P--P+ T 2 P- 



P++P- 



P+ + P- 



dS 



2( P+ + p_)J St * 



+ 



P+ + 



— [ V„t (-A 
P- Js t ~ 



(-A St tf)^(-A St )-*vi 
s t rt)^(-As t )hl 



dS 



dS 



then integrate by parts in these last two integrals estimating Dv± in L°°(St) and the remaining |fc — \ derivatives on v\ in L 2 {S t ). 
For the terms involving the second fundamental form, (A. 9) gives 

\A St (Il ± (vl,vl)-V 2 k ± (vl,vl)\ H s k _ l(St) <Q. 
Since V 2 k± = V„t V^t k± — T> v t v± ■ Vk± and the last term in this sum is lower order, we get 



P++P- 



-V„tV„tk„ 



P+ 



P++P- u+ u+ 



dS 



<Q. 



Using the same previous argument we can commute one of the factors V^t and move it outside to obtain 



dt 



E x + 



P- 



P+ + P- 

P- 
P+ + P- 



[ X7 vT (-A St viW(-A St Af) k - 2 V v T K _dS 
J St 

■ f v vl {-A St vi)rt{-A St rt) k ~ 2 v vl K + ds 

J St 



(B.9) 
(B.10) 



/ vi(-A St M) k K+ dS 
J s t 



< Q. 



Now, thanks to identity (A.7) 



\-A St vi--D t+K+ \ H 3 k _ 2{St) <Q 
so that we can substitute D t+ «+ to -As t u| in (B.9) and (B.10). The usual commutator estimates imply 

d 



dt 



V v t D t+ K + JV(-A St Af) k 2 V vI K_dS 



P+ + P- Js 
P+ 



P+ + P- J s t 



V v rT) t+K+ Af^A St Af) k \ v rK + dS 



< Q 



and (B.7) follows. 

Evolution of E 2 : As before commutator estimates (A. 13) and (A. 14) give 



±E 2 - f K + J\T(-As t Sf) k 'Dt+K+dS 
at J s t 



< Q 
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and in view of (A.7) and (A. 10) we obtain 

d „ 



dt 



n + N(-A St M) k \-A St )v^dS 



St 



< Q 



K+ Af(~A St M) k ^^K+dS 



The same commutation argument previously adopted shows that 



f n + M(-A St M) k \ vl K + dS-\f V„x A^-A^A/vV 

J St J St 



dS 



< Q. 



Integrating by parts and estimating DvJ_ in L°°(S t ) and the remaining |fc — 1 derivatives on k + in L 2 shows that this last integral is 
bounded by Q. Finally use the self-adjointness of M and As, to obtain 



d F 



s, 



v^(-A St SJ) V dS 



< Q. 



(B.ll) 



Evolution of the vorticity oj = Dv — (Dv)*: Commuting repeatedly D t with D and using the identity 



D t w = DD t v - (Dv) 2 - (DT> t v)* + ((Dv)*) = ((Dv)*) - (Dv) 2 = -ljDv - (Dv)\ 



we have 



dt J-j^n 



IDs' 



dx 



-St 



D2 k ~ 1 uj\ dx < C\v(t, 



\u(t, 



l JT3* _1 (R»v.S t ) 



< Q. 



Summing up (B.7), (B.ll) and (B.12) we get the desired cancellations giving (B.5) □ 



B.3 The Energy Inequality 

Integrating in time (B.5) gives 

E(t) - E(0) - E ex (t) + E ex (0) < f Q(\v(s,-)\ Him ^ Ss) , K^Oy*-! ) ds 

J 

for any < t < min{io, t 2 }. We can estimate the extra energy term (B.6) by 



(B.12) 



\E a (t)\ < 



tf*(A St AT)* 1 V„t« + 



< ck ■ v K+ | ff 3 fc _ i(5t) < C|«(t, -)l ff |,-| (R „ xSt) -)l H | fc -| (5t) 

where the positive constant C depends only on the set A . Interpolating v between iji fe ~i and Hi k and k between H^ k ~i and 
H2 k ~ 1 yields 

\^\<Ie + Ci (i + M^- J(RBsSt) ) 

for some integer to where the constant Ci, which includes |/c| 3 fc _ 5 , depends ultimately only on Do and Ao in view of (B.l). Using 
Euler equations (1.16) and lemma 4.1 to estimate the pressure, we have 



p 



1 I v p+IhI-I ( o+) + ^I v p- 



We can then use the Lagrangian coordinate map to estimate 



l^-Cf*-§( R n v S t ) - KMI£$*-J(R»^„) 



< f Q(s)ds 
Jo 



and obtain 



\E ex \<^E + C 1 (l + K0,-)|™3 fc - l(RnxSo) )+^ Q(s)ds< 1 -E + C 2 + Q(s)ds 
where C 2 is determined by D , the set A and I v(0, •) I 3 fc _ 3 Inserting this last inequality in (B.12) we finally obtain 

E(S t ,v(t,-)) <3E(S ,v(0,-)) + C 2 + [ Q(s)ds 

Jo 

for some C 2 as above. Taking \i in (3.1) large enough compared to the initial data concludes the proof of theorem 2.4 ■ 
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